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Electrical response of two diffusive metals is studied when they are linked by a magnetic insulator hosting 
topologically stable (superfluid) spin current. We discuss how charge currents in the metals induce a spin 
supercurrent state, which in turn generates a magnetoresistance that depends on the topology of the electrical 
circuit. This magnetoresistance relies on phase coherence over the entire magnet and gives a direct evidence for 
spin superfluidity. We show that driving the magnet with an ac current allows coherent spin transport even in 
the presence of U(l)-breaking magnetic anisotropy that can preclude dc superfluid transport. Spin transmission 
in the ac regime shows a series of resonance peaks as a function of frequency. The peak locations, heights 
and widths can be used to extract static interfacial properties, e.g., the spin-mixing conductance and effective 
spin Hall angle, and to probe dynamic properties such as the spin-wave dispersion. Ac transport may provide a 
simpler route to realizing nonequilbrium coherent spin transport and a useful way to characterize the magnetic 
system, serving as a precursor to the realization of dc superfluid spin transport. 

PACS numbers: 72.25.Mk, 75.76.+j, 75.70,-i, 85.75.-d 


Introduction .—Understanding spin transport via collective 
magnetic excitations is currently gaining attention m. An ex¬ 
citing frontier explores how analogs of conventional superflu¬ 
idity, as observed in liquid "^He, can be obtained in magnetic 
systems, and how dissipationless spin currents can be realized 
and detected in such systems ||2l. Conventional superfluid¬ 
ity is characterized by a rigid U(l) order parameter, a single 
quantum-mechanical wave function describing a macroscopic 
number of constituent particles. Dissipationless current, be¬ 
ing proportional to the gradient of the U(l) phase, appears 
only in the phase-coherent state with broken gauge invariance. 
A certain class of magnetic insulators, such as ferromagnetic 
insulators with easy-plane magnetic anisotropy, are also char¬ 
acterized by a U(l) order parameter, whose magnitude and 
phase characterize the magnetic order within the easy plane. 
Here, dissipationless spin current (polarized out of plane) is 
triggered by a collective reorientation, i.e., a planar spiraling 
texture of the magnetic order. Theoretical proposals for real¬ 
izing and detecting such superfluid spin transport have been 
put forth for (ferro- and antiferro-) magnetic insulators HSU 
and multiferroic materials ||6l. 

Superfluid spin transport may be detected using a two- 
terminal spin transport setup, in which two spin-orbit coupled 
metals are attached on opposite ends of the magnetic insulator 
and direct (inverse) spin Hall effect facilitates spin injection 
(detection) Ena. While the two-terminal spin conductance 
so obtained can reveal a smoking-gun signature of spin su¬ 
perfluidity—with Gilbert damping spin supercurrent decays 
algebraically over space while spin current carried by inco¬ 
herent thermal magnons decays exponentially over the spin 
diffusion length—such probe requires a study of spin trans¬ 
mission through various sample sizes. Additional signatures 
of superfluid transport obtainable from a single sample is de¬ 
sirable. On another note, a source of difficulty in realizing dc 
superfluid spin transport are magnetic anisotropies that break 
U(l) symmetry crucial for genuine superfluidity. Superfluid 
spin transport in the dc regime requires the U(l) magnetic or¬ 


der parameter to make full 2n rotations within the U(l) plane. 
This dynamics can be quenched by such anisotropies which 
tend to pin the texture. While this pinning may be overcome 
by injecting large enough spin currents |l2|, achieving such 
large spin currents does pose experimental challenges a. 

In this Letter, we theoretically propose measurements based 
on the two-terminal setup that can help overcome these chal¬ 
lenges and thus expedite the realization and verification of su¬ 
perfluid spin transport. First, we show that a nonequilibrium 
spin superfluid state induces unique nonlocal magnetoresis¬ 
tance signatures in the metallic contacts used for spin injection 
and detection. Due to phase coherence over the entire mag¬ 
net, the nonequilibrium superfluid state is self-consistently 
determined by boundary conditions defined at the injection 
and detection interfaces. In the presence of magnetoelectric 
coupling at these interfaces, the charge currents in the metal¬ 
lic contacts define a particular superfluid state, which in turn 
governs a nonlocal magnetoresistance in the circuit. This phe¬ 
nomenon, which we refer to as spin-superfluid magnetoresis¬ 
tance (SSMR), depends intimately on the topology of the ex¬ 
ternal electrical circuit connecting these contacts. The obser¬ 
vation of SSMR would constitute a strong signature of spin 
superfluidity even without sample-size dependence studies. 

In the latter half of the work, we consider spin transport by 
coherent magnetic textures in the ac regime and in the pres¬ 
ence of U(l)-breaking inplane magnetic anisotropy. Even if 
the dc superfluid transport is quenched by pinning, coherent 
ac spin transport can still occur via oscillations of the local 
magnetic order about the pinning potential minimum (without 
the need for full 2n rotations as in the dc regime), thus cir¬ 
cumventing the pinning problem. Spin transmission through 
a magnet of length L shows resonances as a function of ac 
frequency w. For w » wq, being the gap in the spin- 
wave spectrum associated with the inplane anisotropy, the res¬ 
onance peaks occur in intervals of nvfL, allowing one to ex¬ 
tract the spin-wave velocity v. Spin transmission is exponen¬ 
tially suppressed for w « loq, and the position of the first peak 
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FIG. 1. (color online) Schematics of the series (a) and parallel (b) 
configurations, as detailed in the text. 

(at a> - a>o) gives a direct measure of the gap. Furthermore, 
measurements of peak heights and widths allow one to extract 
the effective spin Hall angle and spin-mixing conductance at 
the interfaces. We show that finite-frequency spin waves can 
support spin transmission decaying algebraically as a function 
of sample length, as in spin transmission via dc spin superflu¬ 
idity. Therefore, the ac transport studies should not only serve 
as a simpler route to realizing nonequilbrium coherent spin 
transport, but also as a useful way to characterize the magnetic 
system and a meaningful precursor to the ultimate realization 
of the dc superfluid spin transport. 

General considerations .—Consider a magnetic insulator 
sandwiched by two normal metals as sketched in Fig.[T] Our 
focus is on magnets well below the magnetic ordering temper¬ 
ature, which, in the long-wavelength limit, is characterized by 
a slow continuum variable, the U(l) order parameter n(x, t), 
encoding magnetic state in the (easy) xy plane. Specifically, 
recall that for a ferromagnetic insulator, n corresponds to the 
direction of the local spin density, and the U(l) easy plane 
can generally be defined by the shape anisotropy Q. For 
an isotropic antiferromagnet, n is the direction of the local 
Neel order, and the U(l) plane is defined to be normal to a 
uniform external magnetic field lH. For the axially symmet¬ 
ric magnetic state, the spin density polarized along the z axis 
is a hydrodynamic quantity that is approximately conserved. 
(Its relaxation is in practice governed by spin-orbit impurities, 
which microscopically break the symmetry.) For simplicity, 
we take the normal metals and the interfaces to be identical 
on the two sides. The metals, treated here as diffusive films 
of thickness lying parallel to the yz plane, possess strong 
spin-orbit coupling with an effective spin Hall angle 6 at the 
interfaces. 

The left (/) and right (r) interfaces, located at x = 0 and 
X — L, respectively, break translational symmetry along the x 
direction, but full translational symmetry is assumed parallel 
to the interface (yz) plane. The entire heterostructure can thus 
be described using three coarse-grained hydrodynamic fields; 


the U(l) phase (p(x, t) and out-of-plane spin density s^(x, t) in 
the magnet and the 2D charge current densities in the left and 
right normal metals, j,(f) = (j], /,) and jr(t) = (jl, fi). For 
concreteness, we hereafter focus on an easy-plane ferromag- 
net 0 : n = s/s, in this case, being the direction of the local 
spin density s ^ (s cos ip, s sin tp, s^), where s is the magnitude 
of the equilibrium spin density associated with the magnetic 
order. The formalism is sufficiently general that it can be read¬ 
ily extended to other magnets supporting spin superfluidity; 
it is straightforward to show, in particular, that the case of a 
Heisenberg antiferromagnet is closely analogous 14]. The dy¬ 
namics of an isolated easy-plane ferromagnet is given by 0 

K A 

tp - —+ ah^ , - —tp" - atp , (1) 

s s 

where A and K parameterize the exchange stiffness and the 
easy-plane magnetic anisotropy, respectively, and a is the 
Gilbert damping parameter. The primes (dots) denote differ¬ 
entiation with respect to x (time). Recognizing the second 
equation in Eq. Q a the continuity equation for s^ = sn^, the 
z-polarized spin current (hereafter referred to as simply spin 
current) reads f(x,t) - -Aip'(x,t). 

In the presence of an external electric field E, a uniform 
current-carrying state of an isolated metal is governed by 
Ohm’s law pj(t) - E(t), where p is its (2D) resistivity. In 
the presence of spin-orbit coupling at metal|magnet interfaces, 
current in the metal can induce a torque t on the adjacent 
ferromagnetic moments, and, inversely, the ferromagnetic dy¬ 
namics would induce an electromotive force in the adjacent 
metal. According to spin Hall phenomenology 0 , the torques 
at the left and right interfaces can be written as 

T^r = ±(?7 + X ji^r) X n/,r, (2) 

the upper (lower) sign corresponding to the left (right) in¬ 
terface, and constants 77 and § quantifying the field-like and 
damping-like torques, respectively. Here, n/(f) = n(x - 0, t) 
and ririt) = n{x - L, t). The coefficient for the damping-like 
torque can be related to the effective interfacial spin Hall an¬ 
gle 6 via d = li tan 01 letN 0. By the Onsager reciprocity, 
the torque in Eq. 0 gives rise to an electromotive force ei^r 
in the adjacent metals, thereby modifying the Ohm’s law to 
pjir = Ei^r + £i.r, where 

£lr = ±[(?7 + dniry-)nir] XX . ( 3 ) 

In the following, we will retain only the y components of these 
electromotive forces, as the z components are counteracted by 
an electrostatic buildup along the z axis (supposing the mag¬ 
netic dynamics are slow compared to the relevant RC time of 
the metallic terminals). 

In addition, a physical contact to the adjacent metals gives 
rise to an interfacial contribution to Gilbert damping for the 
ferromagnet. This damping modifies the torques to T/,r ^ 
= Ti^r - X fi;_r, where and is the 

effective (interfacial) spin-mixing conductance. The torques 
tI^ reflect the spin currents entering the ferromagnet at each of 
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the interfaces, providing the boundary conditions for the mag¬ 
netic dynamics in the bulk. The electromotive forces in Eq. (|^ 
(which enter in the modified Ohm’s law) quantify the feed¬ 
back applied by the magnet on the external electric circuit. 
Together these ingredients constitute self-consistent magneto¬ 
electric dynamics, which we systematically address below. 

Nonlocal magnetoresistance .—Our metallic contacts are 
integrated into an external electrical circuit such that in the 
series configuration [see Fig. [TJa)] the two metals have cur¬ 
rents running in the opposite directions, while in the parallel 
configuration [see Fig. [^b)] they run in the same direction. 
For a time-independent E, the hydrodynamic variables take, 
according to Eqs. Q’ a steady-state form: (p{x, t) - f(x) + Q.t 
and rij = const 0, where (,0 s Q is the uniform global pre¬ 
cession frequency of the magnetic texture [to be determined 
self-consistently from Eqs. Q and the boundary conditions 
(see below)]. Matching the torques 0 with the spin currents 
in the magnet near the two boundaries, we ari'ive at the fol¬ 
lowing boundary conditions for f{x): 

-Af(0) ^z-t; , -Af'iD = -i ■ 7 -;. (4) 


For the series configuration, we have ji - -jr - jy. By 
inserting the steady-state form for (f and into Eqs. Q and 
0- the rotation frequency within linear response becomes 


Q = 


§ 

yll +yal2^ ’ 


(5) 


with ja = asL. This result is analogous to that obtained in 
Ref. 0, but now recast in terms of spin Hall phenomenol¬ 
ogy. In linear response, = ±§Q.y, so that the modified 
Ohm’s law gains an additional magnetoresistive contribution 
(p + Pm)j - E, where we obtain \Ei\ - |FJ,.| = E under the as¬ 
sumption of identical metals and symmetrical interfaces, with 




yB + jal'l 


( 6 ) 


The spin superfluidity thus reduces the effective resistivity 
of the circuit, implying that magnetic dynamics reduces net 
power dissipation, for a fixed current. 

For the parallel configuration, ji - jiy and jV = jry, and 
the rotation frequency becomes, within linear response. 




jl - jr 

7^^ + Jal^ 2 


(7) 


Since the electrical circuit is parallel, Ei - Ey Solving for 
the charge currents flowing in the metals, we obtain ji - jr - 
Ejp, so that the magnetic texture is static, Q = 0, a result 
consistent with the mitTor symmetry about the x - Ljl plane. 
The resistivity of the electric circuit is not modified in this 
configuration. SSMR can be distinguished from the spin Hall 
magnetoresistance (SMR) recently discussed in the context of 
ferromagnet-metal interfaces ITTUTSI . SMR generates longi¬ 
tudinal corrections to electrical resistivity of order O^AffltN 
(when tff > An, the electron spin diffusion length) for both 
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FIG. 2. Dynamic transresistance for L = lOOf as a function of ac 
frequency a). Here, we set a = 10“^ and = 10“^. The solid 

line corresponds to the case without inplane anisotropy, while for the 
dashed line we use k/K = Ax lO^”*. 


circuit configurations ||2l[l2l^ contrasting with SSMR, which 
is nonzero only for the series configuration. 

Coherent ac spin transport .—The ac regime can be ex¬ 
plored by driving the metallic contacts by oscillating electric 
fields, i.e., Eiy{t) = Ef^e^'‘^‘. In this section, we introduce 
the in-plane [U(l)-breaking] magnetic anisotropy by adding a 
free energy density term 12 (with k « K), which 

augments Eq. Q to 

K A ,, K 

if - —n^ + ah^, h^ - —tp - aip - ip . (8) 

s s s 

In the steady state, within linear response, the relevant hydro- 
dynamic variables should oscillate at the ac frequency, such 
that ip{x, t) - f(x)e~“^‘ and n^(x, t) - g{x)e~“^‘. The functions 
/ and g can then be obtained using Eqs. 0 and 0. Through¬ 
out this section, we set ji{t) = j{t)y and jr = 0. Assuming 
a large eccentricity of precession, so that the ratio of torque 
coefficients is |i9'/77| » \g/f\, the electromotive force induced 
in the detector (right) contact reads £rit) - id^u)f{L)e^“^\ We 
then evaluate the transresistance as p, = 6^(0/ jit). 

There are two notable lengths scales, which determine the 
loss of spin transmission due to Gilbert damping. Previously, 
it was shown that a spin current cari'ied by the zero-frequency 
mode (the superfluid component) decays algebraically as a 
function of system size L, and that the role of Gilbert damp¬ 
ing becomes negligible for L ^ Lr, = y^^/as 0. For the 
easy-plane ferromagnet, spin current carried by a coherent 
finite-frequency spin wave should decay exponentially at dis¬ 
tances larger than where v = 's/AK/s 

and = K/s are the spin-wave velocity and easy-plane 
anisotropy, respectively. Since it is reasonable to assume 
y^^/s£ 1, we have La « ta. When L » ta, all spin- 

wave modes are strongly damped and the transresistance sig¬ 
nal is exponentially small. We will therefore assume L <si fa¬ 
in the ac regime, a series of resonance peaks appears in the 
(modulus of the) transresistance as a function of frequency 
(see Fig. 0. For zero in-plane anisotropy (i.e., /r = 0) and 
a fixed L, these resonances occur at oj = miv/L, with inte¬ 
ger n, such that the spin-wave velocity can be extracted from 
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FIG. 3. Nonlocal magnetoresistance as a function of system size 
L. Here, we set a = 10“^, !s( = 10“^, and KjK = 4 x 10^“* 

(corresponding to ojoT = 0.02). The solid line is for air = 0.05 > 
WqT, while the dashed line for ojt = 0.019 < WqT. The inset shows 
the semilog version of the same plot. 


the peak intervals (see solid line in Fig. [^. In the presence 
of in-plane anisotropy /c > 0 (see dashed line in Fig. |^, the 
locations of the resonance peaks remain unperturbed for ac 
frequencies well above the gap scale, i.e., oj 3> a>o = 

(see the shaded region in Fig.|^. Spin transmission is expo¬ 
nentially suppressed for at cjq and the first resonance peak 
is shifted to cjq, allowing one to quantify the strength of the 
in-plane anisotropy. The resonant peaks have the square-root 
Lorentzian form, where the height \p,\ and the full-width half¬ 
maximum 6u) (for n > 1) are, respectively, given by 


IP/I 


2yTi + yJ2 ’ 


6(jj - 


2 V3 

stL 


(2r'U|). 


( 9 ) 


The ratio of peak heights between the n > 1 resonances and 
the n = 0 resonance is given by 


ip/i«>i _ 2 


Ip/l«=0 27TJ- -H yal' 2 . ' 


( 10 ) 


The overall factor of 2 and the factor of 2 in the denominator 
arise as a result of finite-frequency spin wave modes having 
antinodes at the two interfaces; this leads to an enhancement 
of interfacial damping (i.e., the factor of 2 in the denominator) 
and an overall increase in the transresistance signal (i.e., the 
overall factor of 2). The spin-mixing conductance, the spin 
Hall angle and the Gilbert damping parameter can all be ex¬ 
perimentally determined for a given device via Eqs. @ and 
( [T0| i. The measured values of dcj and the ratio determine the 
spin-mixing conductance and the Gilbert damping parameter. 
The spin Hall angle can finally be obtained from the measured 
peak height. 

The dependence of transresistance \p,\ on magnet size L re¬ 
flects magnetic losses in the bulk. When the ac frequency 
exceeds the gap (see solid line in Fig. [^, resonance peaks 
are observed up to the scale L ~ but beyond this length 
scale, they become exponentially suppressed (as shown by the 
shaded region in Fig. inset). For L « (a, the peak heights 


decay algebraically, as in the case of spin transmission via dc 
spin superfluidity [cf. Eq. (|^]. When co < coq, spin-wave 
modes cannot be resonantly excited and the transresistance 
decays exponentially with increasing L in a monotonic fash¬ 
ion (see dashed line in Eig.|^. In the exponentially decaying 
region (when either L » or w < oiq), the inverse decay 
length (i.e., the slope of the semilog plot in the Eig. [^in¬ 
set) is given by 


r' 


^ H- (awT)^j*^"^ X 


[cos I , 


CO < coo , 
CO > coo , 


( 11 ) 


where A s (co^ - cOq)t^ and 0 s tan“*(awT/|A|). 

Conclusions .—We studied the phase-coherent dc and ac 
spin transport through a magnetic insulator, sandwiched by 
two strongly spin-orbit coupled metals used for spin injection 
and detection. In the dc regime, we consider spin supercur¬ 
rents carried by both dynamic and static magnetic textures, 
leading to distinct nonlocal magnetoresistances in the met¬ 
als and a direct evidence for the spin superfluid state. We 
also address coherent spin transport in the ac regime, which 
allows one to extract static interfacial properties and the dy¬ 
namic properties of the magnetic bulk, and to realize coherent 
spin transport even when the dc superfluid state is quenched 
by U(l)-breaking anisotropies. 
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